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1. Introduction 

These conference proceedings report on preliminary results from CDF [|l|] and DO [§] investi- 
gating whether the present data allow or preclude the production of hypothetical new quarks which 
decay to final states with a high transverse momentum (pr ) lepton, large missing transverse energy 
(fir), and multiple hadronic jets, having large total transverse energy (Hj), and thus mimicking 
top quark pair event signatures in the lepton+jets decay channel. Previous versions of this search 
were conducted at CDF at lower integrated luminosities [Q, ||]. 

We refer to the hypothetical new quark as t' for brevity, although such a signature could be a 
standard fourth-generation up-type heavy quark in which the splitting between the f'mass and the 
b' mass is less than the mass of the W boson (so that the decay is predominantly to Wq), as well 
as any up-type quark. For the purpose of setting limits we assume that the new heavy quarks are 
produced strongly, are heavier than the top quark and decay promptly to Wq final states. 

Due to the variety of theoretical models predicting similar signatures as well as the number of 
free parameters within each model, an a priori method was established to estimate the significance 
of a potential excess of events without attributing the excess to a particular new physics model. In 
our analysis no significant excess is observed; we thus set a limit on the fourth-generation t' quark 
pair production cross section (times branching ratio of t' — > Wq). 

2. Theoretical Motivation 

There are a number of new physics models predicting heavy quarks with masses above the one 
of the top quark whose decays produce event signatures similar to those from top quark decays. 
The interested reader is kindly asked to refer to the full conference notes and reference therein for 
further details and list of references [jl], Here we only outline a few cases of interest to motivate 
this search. 

One of the simplest extensions to the standard model with three generations is a fourth chiral 
generation of massive fermions. Although not popular historically, such an extension is in good 
agreement with electroweak precision data. To avoid the Z— > vv constraint from LEP I a fourth 
generation of neutrinos must be heavier than 1/2 the mass of the Z-boson. Similarly, to avoid the 
LEP II bounds, a fourth generation of charged leptons must have m > 101 GeV. At the other end of 
the spectrum, sizable radiative corrections means masses of fourth generation fermions cannot be 
much higher than the current lower bounds; the masses of the new heavy quarks t' and b' should 
thus be in the range of a few hundred GeV, within the reach of the Tevatron collider. In most models 
a small mass splitting between the t' and b' is preferred such that the t' decays predominantly to 
Wq where the q is one of the known standard model down-type quarks. The present bounds on the 
Higgs are relaxed in the presence of a four-generation model leading to Higgs masses as large as 
500 GeV. Furthermore CP violation is significantly enhanced to a magnitude that might account 
for the baryon asymmetry of the universe. 

3. Data Samples, Monte Carlo Simulation 

The data are collected at the CDF and DO detectors at the Fermilab Tevatron proton-antiproton 
collider at a centre-of-mass energy of = 1.96 TeV The CDF and DO detectors are described 
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elsewhere The latest analysis from CDF (DO) uses 4.6 fb _1 (4.3 fb _1 ) of integrated lu- 

minosity. The events must have fired one of several trigger conditions, all but one requiring an 
isolated electron or muon with a high transverse momentum. CDF increased the signal acceptance 
by incorporating muons collected with a f.j + jets trigger. 

The t't' signal is generated with PYTHIA 6 [0] (v6.216 for CDF and v6.409 for DO). For the 
CDF analysis the branching ratios are not modified, for the DO analysis the branching ratio to Wb 
is set to 100%. In both cases, as no b-jet identification is applied, the results are applicable to 
both heavy quarks decaying to a W and light or a heavy quark. We assume that the t' is a narrow 
resonance and as such the intrinsic width will be much smaller than the detector resolution and 
therefore the exact value of the resonance width does not affect the analysis. Various mass points 
were generated for both analyses in order to cover the t' mass range between 200 and 500 GeV. 

The background consists mostly of 

• tt which is modeled with ALPGEN [J8j] and use PYTHIA for the parton shower for DO and 
PYTHIA 6.216 for CDF, 

• W + jets which is modeled with ALPGEN (v2. 10' for CDF and v2. 1 1 for DO) and use PYTHIA 
(v6.325) for the parton shower with a jet-matching algorithm following the MLM prescrip- 
tion [g], 

• Multijet events are those events were one of the jets fakes a lepton; this contribution is mod- 
eled using data. The methods used by CDF and DO to obtain the rate and shapes of the QCD 
background are described in the public notes [jl|] [^J. 

Smaller contributions from dibosons, Z + jets and single top are also taken into account. 
4. Event Selection 

While the specifics of the cuts differ between the two experiments, the principles of the cuts are 
similar. In addition to the trigger selection mentioned in the previous selection, events are selected 
by requiring a single isolated high-/? r electron or muon, high missing transverse energy, and four 
or more jets in the event (with slightly higher cuts on the leading one or two jets). 

In the CDF analysis two additional classes of cuts are applied. The first set of cuts are on the 
transverse mass of the W and the missing transverse energy significance to reduce the contribution 
from QCD events which are notoriously difficult to measure, especially the tails where the signal 
is expected to be important. Such cute are also implemented in the DO analysis. The second set 
of cuts are targeted at cleaning up the tails of the distributions and to as much as possible remove 
mis-measured objects causing large values of missing transverse energy. These are cuts that are 
applied to the high-/?r electrons, muons and jets and ensure that the missing energy is not back 
to back with the lepton and not roughly collinear with the jets. These two sets of additional cuts 
increased our confidence in the modeling of the tails of the distributions used to fit the t' signal. 
The additional clean-up cuts used by CDF are described in the latest public note [|l|]. 
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5. Variable Selection and Fitting 

We use two variables to distinguish the t' signal from the standard model backgrounds: the 
reconstructed mass of the top it') for CDF (DO) (M rec0 ) and the scalar sum of the transverse mo- 
mentum of all jets, the charged lepton and the missing transverse energy {Hj). Both variables 
exploit the large mass of the t' quark. 

For both experiments the mass is reconstructed using the minimization of a % 2 tit describing 
how compatible the event is with a top (?') decay hypothesis. While the exact form of the % 2 
functions differ, the concept and performance is similar between the two experiments. 

The CDF analysis further separates the data into four different sets of distributions: events 
with exactly 4 jets and those with 5 or more. Each is further separated into events where the x 2 
value of the best fit is less than or greater than 8. The fit is performed simultaneously to all four 
sets of distributions. 

The two-dimensional binned distributions of Hj vs M reco is used to test for the presence of a 
t'f signal in the data. As no signal was observed, the 95% confidence level upper limits on the 
production cross section are computed. 

The CDF analysis calculates a likelihood as a function of the signal production cross section 
and uses B ayes' Theorem to convert it into a posterior density in a t 'p. This is then used to set an 
upper limit or measure the production rate of the signal. The production rate for W+jets in each 
jet bin are two free unconstrained independent parameters of the fit. Other parameters such as the 
tt production cross section, lepton ID data/MC scale factors and integrated luminosity are related 
to systematic errors and treated in the likelihood as nuisance parameters constrained within their 
expected (normal) distributions. We adopt the profiling method for dealing with these parameters, 
i.e. the likelihood is maximised with respect to the nuisance parameters. Systematic uncertainties, 
such as the jet energy scale, the W+jets Q 2 scale uncertainty of the initial and final state radiation 
uncertainty enter the likelihood as gaussian constraint penalty terms. Outside the interval [-1,1] the 
penalty terms are extrapolated as a linear function of the jet energy scale parameter. This treatment 
is called vertical template morphing. 

The DO analysis uses as a test statistic the likelihood ratio L = —llog^p^-, where Pg (Ps+b) 
is the Poisson likelihood to observe the data under the background only (signal plus background) 
hypothesis. For the background only hypothesis, three components are fit to the data: tt production 
constrained to its cross section, W+jets constrained to the number of expected events with their 
uncertainty, all other backgrounds added together in their expected proportions and constrained 
with an overall normalization that floats freely. To determine the cross section limit, the CL$ 



method is used JlOfl- The systematic uncertainties are treated in a similar way as for CDF. 

The CDF analysis uses 28 bins for Hj and 18 for M reco with the overflow bins defined as events 
with Hj > 800 and M recu > 500 GeV. Due to limited MC statistics, an algorithm was developed 
to merge contiguous bins with low total MC statistics to form super-bins used in the likelihood. 
This procedure by construction deteriorates the sensitivity of the new physics signal, however it 
eliminates the abnormalities from sources such as due to bins with zero predictions and thus provide 
reliable observed limits. The exact algorithm and criteria used to define the merging are described 
elsewhere [fiTl. 
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6. Results 

The sensitivity of the methods are tested by drawing pseudo-experiments from standard model 
distributions. The range of expected 95% CL upper limits within one and two standard deviations 
are shown in Fig. [I] as the colored bands. The straight line falling curves represent the theoretical 
predictions. The left plot shows the results from CDF, the right one is from DO. From these figured 
it follows that given no t' presence, the methods are on average sensitive to setting upper limit at t' 
masses of 372 GeV (CDF) and 330 GeV (DO). We perform the analysis fit on data and determine 
the upper limits on the t' signal. The red curves show the final results. 

Based on these results we conclude that there is no significant excess of observed events in 
either experiment and we can exclude a t' mass below 335 GeV at CDF and 296 GeV at DO, given 
the true top mass is 172.5 GeV. 



DO Run II preliminary 4.3 fb 1 




t' mass (GeV/c 2 ) r mass [Gev] 

Figure 1: Observed and expected cross section limits on t' —> Wq compared to theory for CDF (left) and DO 

(right). 
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